The Leading Off-Diagonal Correction to the Form Factor of Large Graphs 
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Using periodic-orbit theory beyond the diagonal approximation we investigate the form factor, 
K(t), of a generic quantum graph with mixing classical dynamics and time-reversal symmetry. We 
calculate the contribution from pairs of self-intersecting orbits that differ from each other only in 
the orientation of a single loop. In the limit of large graphs, these pairs produce a contribution 
— 2t 2 to the form factor which agrees with random-matrix theory. 
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One of the fundamental open questions in the quantum 
theory of classically chaotic systems is whether the fluc- 
tuations in the spectra follow the predictions of random- 
matrix theory (RMT) ll| and if so, precisely under which 
conditions. An interesting spectral statistics to consider 
is the form factor K(t), which is the Fourier transform 
of the spectral two-point correlation function and de- 
pends on a dimensionless time r measured in units of 
the Heisenberg time. We will concentrate on systems 
with time-reversal symmetry, the appropriate ensemble 
of random matrices is the Gaussian orthogonal ensemble 
(GOE) and the corresponding form factor is B 




K GO e(t) = 2r-rlog(H-2r) (0 < r < 1) 

= 2t-2t 2 + 0(t 3 ). 



(1) 



There is ample numerical evidence JG[ for the conjecture 
that in the semiclassical limit, Eq. |l]) correctly describes 
the spectral fluctuations of any "typical" system with 
chaotic classical limit, but there is no complete analyti- 
cal theory. A natural starting point for such a theory is 
the Gutzwiller trace formula || expressing the spectral 
density in terms of classical periodic orbits (POs). How- 
ever, until recently, the only method for dealing with the 
resulting double sum over POs was Berry's diagonal ap- 
proximation [0 which rests on the assumption that clas- 
sical POs are uncorrelated unless related by exact sym- 
metries. This can be justified only in the limit t — > 0, 
and the diagonal approximation results in the leading- 
order term It but fails to reproduce any higher-order 
corrections in Eq. (fy). 

The analogy to diagrammatic perturbation theory for 
disordered systems suggests that the leading correction 
to the diagonal approximation is related to correlations 
between pairs of orbits which differ only in one self- 
intersection || g] (see Fig. [j]). These correlations were 
recently calculated for certain Ricmannian surfaces with 
constant negative curvature M, p|, and they were found 
to reproduce the — 2r 2 term in Eq. (|]). 

In this letter, we study the contribution from self- 
intersecting orbits in a very different class of systems: 



FIG. 1: In (i) we show two loops on a graph intersecting 
at a vertex a. The dashed arcs indicate trajectories between 
vertices a, b and c, d, respectively. They may contain interme- 
diate vertices which are not necessarily pairwise different. In 
(ii) we show all orbits which traverse each of the loops exactly 
once. Dashed arrows between orbits indicate pairs that were 
counted in the diagonal approximation. Solid arrows indicate 
pairs missed by the diagonal approximation but counted by 
our method. 



quantum graphs B. In these systems two-point corre- 
lation functions can be expressed as combinatorial sums 
and therefore allow a treatment with fundamentally dif- 
ferent analytical tools J10|, O, Q. Our result can be 
summarized as follows. The contributions to the form 
factor from self-intersecting orbit pairs differing in the 
orientation of a single loop are calculated analytically. 
We find that they add to — 2t 2 if correlations in the clas- 
sical dynamics of the graph decay sufficiently fast (Eq. (Q) 
below) . The role of the unitarity of the quantum dynam- 
ics is particularly transparent in our derivation and clas- 
sical information enters only via the Frobenius-Perron 
operator — a concept which is not specific to system or 
representation. However, the most important aspect in 
which our result extends Refs. [0, |) is that it applies to 
generic graphs and not only to uniformly hyperbolic sys- 
tems where all POs share the same Lyapunov exponent. 

We consider a graph with N vertices and B directed 
bonds. The quantum dynamics on the graph is defined 
by the energy dependent unitary bond-scattering matrix 



S, which can be interpreted as the discrete time-evolution 
operator on the B dimensional Hilbert space of directed 
bonds. The matrix element S m >i t i m = o„ l , m e lfeLmi pre- 
scribes a transition amplitude from bond to — > I to bond 



I 



This amplitude consists of a phase kL m i, where 



L m i is the length of the bond m —> I, accumulated dur- 
ing the free propagation with wave number k from ver- 
tex to to vertex I, and an element of the unitary vertex- 
scattering matrix a^ describing the vertex I g|. Matrix 
elements which do not correspond to any bond-to-bond 
transitions are zero. Solely for the purpose of a transpar- 
ent presentation we assume that for each pair of vertices 
to, I there is at most one directed bond to — v I. To ensure 
the invariance of the system under time reversal we re- 
quire the matrices a"' to be symmetric and L m i = Li m 
for all pairs of connected vertices l,m. Otherwise differ- 
ent bond lengths are assumed rationally independent. 
The form factor of a graph can be defined as p0[ 



K<< B \T) = B- 1 (\trS t (k)\'' 



(2) 



where r = t/B and the average is over the wave number 
(•)fe = lim K _>oo kT 1 L dfc(-). We will be interested in the 
limiting behavior K(t) = liniB^oo K (t), i- e. we first 
fix r and consider larger and larger graphs, a procedure 
corresponding to the semiclassical limit. Then we look 
at the expansion of K(t) around r = 0. 

Associated with the unitary matrix S is the dou- 
bly stochastic matrix M with M m ti^ rn = \S m ii,im\ 2 = 
l°m'ml 2 - This matrix defines a Markov chain on the 
graph which represents the classical analogue of our 
quantum system || |]|. M can be considered as the 
Frobenius-Perron operator of this discrete classical dy- 
namics. Physically, [M'] m /|» jm i is the classical probabil- 
ity to get from bond (I — > m) to bond (V — > m') in 
t steps. We will assume the Markov chain to be mix- 
ing (and consequently also ergodic) , such that any initial 
distribution approaches for long time an equidistribution 
over all bonds, 



lim M, 



ra'l'ml 



l 

B 



(3) 



for all pairs of bonds m—*l,m'—*l'. Eq. (g) is satisfied if 
the eigenvalue unity of M , corresponding to the invariant 
equidistribution, is the only eigenvalue on the unit circle 
while all others are inside. For a given finite graph this 
is a very weak condition. However, we consider a family 
of graphs for which the limits of large graphs N — > oo 
(and, as a consequence, B — > oo) and large times t — ► oo 
are connected by t/B = r being constant: the size of 
the graph sets the time scale for the decay of classical 
correlations. Therefore we make Eq. (0) more precise 
and require as a sufficient condition for the results to be 
derived below 



lim B 2 N 3 (B) max 



M, 



m'l'ral 



l 
B 



= 0. (4) 



This is only slightly stronger than the recent conjecture 
by Tanner |Q, that RMT behavior results in graphs if 
the spectral gap between unity and the magnitude of the 
next biggest eigenvalue vanishes slower than the inverse 
matrix size 1/B. At the end of this letter we give two ex- 
amples, in which equidistribution is approached at least 
exponentially fast as B — * oo, such that our results do 
apply. 

Expanding the trace of the matrix powers in Eq. (0) 
and performing the average over k, we can represent 
the form factor in terms of pairs of returning trajec- 
tories p, q on the graph which share the same length 
K( B \r) = B- l Y, P , q A p A* q 5 Lp , Lq §, 0. Here p,q de- 
note sequences of t vertices p — [pi, • • • ,p t ], such that for 
any i there is a bond Pi — » Pi+\. The length of the corre- 
sponding trajectory is L p = L P1P2 + L P2P3 + . . . + L PtPl 

and the amplitude is A p = cr p ^] P2 cr Pl * P3 ■ ■ ■ a Pt _ x lPl . If the 
trajectories p, q share the same length, then any pair of 
cyclic permutations of p, q will do so as well. All return- 
ing trajectories p related by cyclic permutations form an 
equivalence class P which represents a periodic orbit. As 
t — > oo, only an exponentially small fraction of the set of 
orbits are multiple repetitions of shorter orbits, ignoring 
these we write the form factor as a double sum over POs 



t 



^ B \r)= b -Y J A P A* Q 5 Lp ^ 



Lq 



(5) 



Due to the rational independence of different bond 
lengths all contributions to Eq. ([|) come from pairs of 
POs which make the same number of transitions between 
any given pair of vertices to, I (i. e. m — > / or I — ► to). 
In particular, the pairs P, P and P, P, where bar denotes 
the time-reversal, contribute to (j5|) for any P since time- 
reversal symmetry implies Lp = L-p- (and Ap = A-p). 
These contributions are collected in the diagonal approx- 
imation to the form factor 



r(B) 



t l 



^dfagM = bH^M 



+ A p A^)^2t, (6) 



which reproduces the linear term in Eq. (|lj). To eval- 
uate the sum (||) we used that according to Eq. (||) 
J2p\ a p\ 2 = t^tvM^^ t~ l . We also neglected self- 
retracing orbits P = P, since they represent an exponen- 
tially small fraction as t — > oo. 

It is the purpose of this letter to go beyond the diag- 
onal approximation and investigate the contribution of 
pairs of orbits which differ from each other at one self- 
intersection. This is illustrated by Fig. ffl: (i) is the exam- 
ple of two loops intersecting at a vertex a. (ii) shows the 
four orbits which traverse each of the two loops once and 
thus have the same length. The arrows between these or- 
bits indicate the possible pairs whose contribution should 
be counted in (g); the dashed arrows correspond to pairs 
already counted in the diagonal approximation and the 
solid arrows to those that we have to sum here. 



(1) 
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FIG. 2: Two conditions make the standard code unique: (i) If 
the vertex separating the two loops is ambiguous, we choose 
it such that the length of the loop with unchanged orientation 
is maximal, e. g. in the top pair of orbits we take a. (ii) The 
loop with unchanged orientation must not be time-reversal 
invariant (self-retracing), e. g. the bottom orbits are excluded. 



The two orbits in any pair that is of interest to us 
differ in the orientation along one of the cycles from 
vertex a to a. To count such pairs we define a "stan- 
dard" symbolic code for the orbits of the pair P — 
[pi,P2, ■ ■ ■ ,Pt] = [a, h, a, h] and Q = [a, h, a, l 2 ], where 
h = [P2, ■ ■ ■ ,Pt'-i] and l 2 = [pt'+i, • • • ,Pt] represent se- 
quences of t' — 2 and t — t' vertices, respectively, and 
h = [pt, ■ ■ ■ iPt'+i] stands for the reversal of the sequence 
l 2 . Alternatively, a pair can also be specified by the code 
of orbit P and the position t' of the second occurrence 
of a. For example, the standard code for the pair of 
top right and bottom right orbits of Fig. |](ii) is P = 
[a,a,...,b,a,d,...,c] and Q — [a,a,...,b,a,c,...,d]. 
We stress that the orbit P will in general contain more 
than one self-intersection. It is the pair P, Q which sin- 
gles out the particular one at a. Nevertheless, the code 
for the pair P, Q is not yet unique. For example, the pair 
shown in Fig. @ (i) has two possible codes: 

P = [a, j3, a, 6, j3, a, c, d], Q = [a, (3, a, 6, f3, a, d, c] 
or P = \J3, a, b, /3, a, c,d,a], Q — [/3, a, b, 0, a, d, c, a] 

with t' — 6 and t' — 4, respectively (remember that 
cyclic permutations of the code are irrelevant). To avoid 
this ambiguity we require (i) pt'+i ^ Pt for the standard 
code which excludes the second of the two representa- 
tions shown and requires c ^ d for the first one to be 
valid. Moreover we assume (ii) l\ 7^ ^1 since otherwise 
Q = P such that the orbit pair was already counted in 
the diagonal approximation. Augmented with the condi- 
tions (i), (ii) the standard code is a unique representation 
of all orbit pairs we have to count. Thus we have the con- 
tribution of such orbits to the form factor 



where the sum is over P = [a, a,p 3 , . . . , b, a, c, . . . , d] (in 
analogy to the notation of Fig. [H we write a, b, c, d for 
P2,Pt'-i,Pt'+i,Pt)- This sum is subject to the restriction 
l\ =/= l\ indicated by a prime and also to all restrictions 
implied by the connectivity of the graph (for some se- 
quence of vertices P there might be no orbit). 

Now we compute the contribution ApAU of an orbit 
pair P, Q. The amplitude Aq differs from Ap only due 
to the vertex-scattering matrix elements picked up at the 
self-intersection point a. P gets a contribution a," ci™ 
from the two encounters of vertex a, while Q has the com- 
bination of adjacent vertices exchanged crca <J bd ■ Thus 



A A* — \A \ 2 \ A |2„.(«)„.(°0_(c0*_(a0 
ApA Q - \A h \ \A h \ a da a bc a' ca > a : ■ 



i«i 



(8) 

where A\ i is equal to the product of matrices a over 
the sequence k. After inserting this result into Eq. (M) 
we can sum over the variables pa,. . . ,pt'-2 and get 



■>P3,—,Pt>-2 

over p t '+2, ■ 
K[ B \r) 



\A 



M 



t'-2 



After a similar summation 



*i I 1 * ± aa,ba 

,Pt-i we obtain 

t'— 4 Oi.a.b.c.d 



X(T da °bc a ca a bd 



*&(T). 



(9) 

Performing an unrestricted summation over p 3 , . . . ,pt>-2 
we lost track of the condition l\ ^ h, therefore the con- 
tribution from such self-retracing orbits is removed by 
subtracting K\ srt (r) which will be calculated below. 

Note that for large t in Eq. (||) either t' — 2 or t — 
t' will be large for any t'. Thus we partition the sum 
over t' into two intervals, Yh~=a = I2t'=4 + X)t'= 2 [t/2]+i> 
where [•] stands for the integer part. We use Eq. (g) 
to argue that in the first interval M*~* = 1/P, while 
in the second interval M l ~ 2 = 1/ B. Calling Ki t s(h) the 
contribution of the interval with Small (Large) t', we have 
Ki = K ltS + K hL - K ltSTt . Substituting M^J a = 1/B 
into i^i.L we evaluate the sum over a and use the unitary 
of a^ to show that K\ t L = 0. There are corrections to 
this result due to the deviation from equidistribution of 
M* ~ 2 . However so long as the stronger mixing condition 
given by Eq. (Q) holds these corrections can be shown to 
vanish as N — ► 00. To show this we use the unitarity 
of a and the fact that all matrix elements of M t ~ t are 
bounded from above by unity. 

1/P, then the sum over c 

(a) (a)* 
. be ° ca 

since the matrix a^ a ' is both unitary and symmetric. 
Hence Ki^ simplifies to 



For K ljS we set M^' 



can be written as J2 P a bc a \a (1 - Scd) = 5 ab -w bd „ da 



K 



(B) 
1,S 




E M LLw ( £\ 2 



t— 2 



K i ^ = rEE^mm 



(T) 



t'=4 P 



E 



M t'-2 I (o),2| («)|2 
lvl aa,ba \"da I \"bd 



,.b,d 



(10) 



(«)|2 



In the first sum inside the brackets we can sum \a^ 
over d to produce one. In the second sum we substitute 



r ( Q )|2 _ 
'da I — 



M d(Xiaa and \a { 

t'-2 
aa.ba 



J2a,b M da,aa.Ml n ^ n M ba , nl , 



'" 2 ~~ Mb aia d and then use 
We arrive at 



ML ad ■ 



[t/2] 



k (b) _ _£_ y> y> (m*'- 2 - m*' \ 

1,S ~~ B 2 Z^ Z^ \ M ab,ba m ab,ba) 

t'=4 a,b 

= ^Y (m 2 + M 3 - M^ 2 !- 1 - Ml*/ 2 !) (11) 

_D ^ — ' V / ab.ba 



Finally we compute the contribution of i^ J t (r) to 
Eq. (|) from orbits with h = h (Fig. |(ii)). Such or- 
bits can be described as follows: To the left of vertex 7 
they must have the same structure as one of the orbits in 
Fig. 0(ii) . To the right of 7 they can take any path which 
comes back to 7 (with the exception of completely self- 
retracing paths which are exponentially few in number 
and can therefore be ignored). The sum over all these 

ft-3 



t-2 



configurations to the right of 7 yields M s s or M 



for even or odd t' , respectively. The contribution of the 
transitions from 7 to the left can also be expressed in 
classical terms, resulting in 



K (B) _ 



7,5 



Ml* -m; 



t-2 



~1o,Sf lvl &i,~fo 



M: 



t-3 

~fS,5~f 1V1 S*y t yS 



M 



(12) 



Matrix elements of large powers of M can again be ap- 
proximated by B~ x both above and in Eq. (|ll|). Then 



we see that the first two terms in Eq. (11) are exactly 
canceled by K[ S1 . t (r) . The remaining terms in Eq. (Ill]) 
give our main result 

K x {t) = -It 1 . (13) 

Eqs. (^), ( J13] ) show that the contributions to the form 
factor of the orbits discussed here precisely reproduce 
the first two terms of the RMT result, Eq. (g). Note, 
however, that we have not discussed whether other cor- 
related orbits could produce contributions of the same or 
lower order. This, we believe, can be shown by consider- 
ing more and more complicated pairs of orbits, eventually 
deriving a complete expansion of the form factor in pow- 
ers of r. 

Finally we briefly discuss two specific models, to which 
Eq. ( pL3|) applies. We consider complete graphs with N 
vertices and B = N 2 directed bonds connecting every 
pair of vertices including a bond from a vertex to itself. 
The first model we consider has "discrete Fourier trans- 
form" (DFT) boundary conditions at each vertex p4[ . 
The second has Neumann boundary conditions H 



r(0 



N -l/2 e 2nimn/N J JJj ( DFT ) 

2N~ l — S mn (Neumann) 



(14) 



The DFT graph is a uniformly hyperbolic model: the 
weights of all classical orbits with the same length are 



equal. We have M l ahcd = 1/N 2 for t > 1, i. e. all 
correlations vanish immediately and Eq. (0) is an iden- 
tity. For the Neumann graph, which is not uniformly 
hyperbolic, the situation is more complicated: Here 
we have M 2 J cd = F T (a,c)F T (b,d)/N 2 and M 2 ^ 1 = 
F T +x(a, d)F T (b, c)/N 2 with F t (a, b) = l + {6 ab N- 1)(1 - 
4/7V)*, but substituting t = N 2 t and using (l-4/N) N -> 
e -4 as N —j- 00 it is easy to see that M* converges for 
large graphs exponentially fast to equidistribution. 
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APPENDIX: DERIVATION OF EQ. (§) 



Applying the same inequality to the sum over b we obtain 



We include this appendix (which is only in the nlin.CD 
version of this paper) in the hope that it will interest the 
reader. In the paragraph below Eq. (g) we noted that if 
the condition in Eq. (|j) is satisfied then K\l = 0. Here 
we show how we arrive at this condition. 

We can write K\ l as 



Ki, 



D 



E E K 



cdoi 



M 



t'-2 
aa,ba 



D 



t' = [t/2\ + l a,a,b,c,d 

x (a) («)*(") (a)* /I _s: \ 
X a bc °bd °da a ca I 1 °cd)- 



To estimate it from above we notice that max 
decreases with t (we take the maximum over a, c, d); and 



M ac,da 



so does max 

1 



mi 



-l/B 



Thus 



lfi f L|<-BVniax|M^ d .| 



M 



[t/2]-l 



1 



aa,ba 



B 



max 

a,a,b,c,d 

Next we apply Cauchy-Schwartz inequality to the sum 



e k 



(15) 



over a 



/ , \"da a ca 



< 






EW 



= i 



(16) 



l-Ki.il < --B 2 r 3 max | M^ cdQ | max 



M 



[i/2]-l 
cm, 6a 



i 
B 



1. 



a,c,d 

(17) 

where the last sum is over all vertices a and over vertices 
c, d adjacent to a. Thus to ensure that |-Ki,i| — > as we 
take the limit B, N, t — ► cxj (r = £/B fixed) we demand 
that 



B 2 max |M 2 cda | max 



71/ 



Bt/2 



B 



E^-*°> ( 18 ) 



where w Q is the valency of vertex a. To make the es- 
timate less dependent on the graph's topology we note 
that M 2 j„ < 1 and «„ < iV, then 



B^rnax|M^V-l/B|^0, 



(19) 



as B (and iV) tend to infinity. Thus we arrive at the 
condition we gave in Eq. (Q). Note a similar calculation 
shows that this is also a sufficient condition for correc- 
tions to Eq. (ft0|) to vanishes. 



